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Abstract - This research reviews the effects of a standard deviation on time response and frequency response of Gaussian filter. In the research, 
standard deviations of 0.5, 1, 1.5, 2, 2.5 and 3 were applied to normalized Gaussian filter in time domain and frequency domain. From the result of the 
simulation, it was observed that in the time domain, the peak amplitude of the filter increased with a decrease in the standard deviation and the 
execution time decreased with a decrease in standard deviation. In the frequency domain, the peak amplitude is constant for all the standard deviation 
and the frequency bandwidth decreased with increase in standard deviation. The simulation was carried out using Matrix Laboratory (MatLab). 


Index Terms - Frequency Domain, Frequency Response, Gaussian Filter, Peak Amplitude, Standard Deviation, Time Response 


1.1. Introduction 

The Gaussian filters find widespread application in signal 
processing because Gaussian masks have non-negative 
values (h (t) 2 0) for the range -- < t < ~, and, therefore, do 
not invert the filtered signal. It is one of the best signals 
filtering method based on the Gaussian distribution 
invented by a German mathematician, Carl Friedric Gauf. 
Gaussian filter is a non-causal filter which has found 
applications in low-level computer vision [1] and pattern 
recognition which includes, image clustering, anomaly 
detection, 
component analysis [2]. The Gaussian filter is a non-zero 


classification, regression and __ principal 


function that theoretically requires an infinite window 
length te[—--, °°] for precise filter implementation. 

The geometric distribution of the mapped time coordinates 
in kernel space of a Gaussian distribution depends on the 
standard deviation (co). The standard deviation of the 
Gaussian function plays an important role in the filter 


e Hyginus udoka Eze is a masters’ degree holder in electronic 
engineering in University of sila Nsukka. Nigeria 
E-mail: ezehyginusudoka@gmail.com 

e Martin C Eze is currently pursuing Phd pee in electronic 
engineering in University of Nigert Nsukka, Nigeria, 
E-mail: martin.eze@unn.edu.n 

e Chidume S Chidiebere is a masters’degree holder in electronic 
engineering in University of Nigeria, Nsukka. Nigeria. 
E-mail:chidumechidieberesunday@gmail.com 

e Bassey Okon Ibokette is a masters’ degree holder in lectronic 
engineering in University of Nigeria, nsukka. Nigeria. 
E-mail;basseyibokette@ gmail.com 

e Micheal Ani is mien pursuing a ae in electronic 
engineering in Univeristy of Nigeria, Nsukka. Nigeria. 
E-mail:micanian4real@gmail.com 

e Uchenna p Anike is a masters’ degree holder in electronic 
engineering in Univeristy of Nigeria, Nsikka.. Nigera. 
E-mail: anikeuc@yahoo.com 


behaviour. The size of a Gaussian filter determines the 
standard deviation needed for effective filtering. The o 
needed is directly proportional to the filter size. The 
selection of an appropriate value of standard deviation (c) 
for different Gaussian filter sizes, is known as tuning [2]. 
For optimum tuning, the value of o must not be too big or 
too small, but must minimize intra-class variation and 
maximize inter-class variation [4] and at the same time 
reduces the execution time. The standard deviations for 
different Gaussian filter sizes are obtained from (1) where 


M is the size of the filter. 
ee ae i 
o= 6 ( ) 


Equation (1) shows that the standard deviation value for 
optimum performance increases with the length of the 
filter. It can be shown that large filter size needs the large 
standard deviation, for effective filtering. The 
comprehensive statistical properties of the Gaussian filter 
create an approach for the Gaussian filter tuning and 
design. 

When the filter size is constant, the variation of standard 
deviation affects the bandwidth, execution time, and 
amplitude of the filter. A filter with a narrow bandwidth is 
used where a very narrow frequency band is needed from a 
broad spectrum. On the other hand, Gaussian filter with 
large amplitude is needed when the input signal is highly 
attenuated. For a situation where short execution time is 
needed, the priority is to choose small standard deviation. 
The execution time of Gaussian filter for 0 < t < ko region 
is calculated usingt, =ko where k is integer which 
determines where the filter is truncated. For -—ko <t < ko 
region, the execution time is given byt, = 2ko. Table 1 
shows how execution time of Gaussian filter is affected by 
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different standard deviation computation expressions for - 
30sts3. o. From the table, it is seen that using 
o= “— technique gives Gaussian filter with better 


performance compared to other expressions. 


1.2. Gaussian Filter in Time Domain 


In the time domain, A-dimensional Gaussian filter, h(t) 
with the mean at the origin is represented as in (2) where t 
is a vector representing time and o represents the standard 
deviation of the filter [3]. 


1 Bre, ti 
h(t,) = (ana) 20 (2) 


When the filter mean is not at the origin, the A-dimensional 
Gaussian filter in (2) is rewritten as shown in (3) with p as 
the mean where the parameter I|t- uw | is called the 
Euclidean distance. 


1 WA i-w? 
h(t,) = (ano) * 207 (3) 


When the value of A in (2) is unity, (2) can be rewritten as 
shown in (4) [1]. 


— 
h(t) = (ano) 2 (4) 


The expression for the error (€(t)) when (4) is truncated 
within —ko < t < ko region is given by (5). 
29 ko 
£6) = ee bom Soyg MCat s) 
f(t) dt 
The 1-dimensional Gaussian kernel is applied in the 
filtering of 1-dimensional signals. 


1.3. Gaussian Filter in Frequency Domain 


The Gaussian function is a special function that occurs 
frequently in mathematics. It is a special function because, 
in the frequency domain, it is Gaussian. The Gaussian filter 
is expressed in the frequency domain using Fourier 
transform. Mathematically, the Fourier transform of a 
Gaussian function is given by (6) where w is the frequency 
(rad/sec). 


t2 
e 207 ,e J dt (6) 


(h(t)) =H = | 
FORO) = HW) = eee | 


The expression in (6) may be expanded as in (7). 


t?2 
e 207 .cos(wt) dt 


Oe ([- 
aa = 


=3 | o3e8 cine) at) (7) 


Since the second integral is odd, the integration over a 
symmetrical range is zero. Therefore, the Fourier transform 
of the Gaussian function is given by (8) [5]. 


wo 


H(w)=e 2 (8) 


The expression in (8) represents a Gaussian filter in the 
w2o2 
frequency domain. From the (8), |H(w)| =e 2— 


and ZH(w) = 0°. 


When a Gaussian filter in frequency domain is truncated 
within —* <w <* the expression for the error (€(@)) is 
given by (9). 

k 
[-H@) dw — f% H(@)dw 


§@) = (9) 


{_. H(w) dw 


1.4. Results and Discussion 


Since the Gaussian function decays rapidly, it is often 
accurate to truncate the filter to enhance the speed of 
operation. In frequency domain, the Gaussian distribution 
is truncated at w = += since the amplitude of Gaussian 


kernel approaches zero when the bandwidth (w) becomes 
more than .. This can be shown in Table2. From Table 2, it 
is seen that atw = +5, the value of H(w) =0.0003 which is 
very close to zero. At any point when w > + 7 the value of 
H(w) is zero. This implies that the filter can be truncated at 
wo=t : with very negligible error. The error in truncating 
Gaussian filter in the frequency domain is calculated using 
the ratio of the normalized area excluded by -* sows * the 
region to the normalized area within 
—c <q <e region. For a normalized 1-dimensional 
Gaussian filter in the frequency domain, the normalized 
area within the —c < w < ~ region is 10.0265. 

The percentage area occupied by a region determines the 
accuracy and the bandwidth of the filter within the region. 
The computation of areas within regions in Gaussian 
distribution is shown in Table 3. In Table 3, the area of a 
Gaussian distribution within the regions, -:%7 < w < »% and 
—= < w <*is 10.0265 each. 

The area occupied by the region -* <a <- accounts for 
68.27% of the total area of the distribution, the region 
= <w <* accounts for the 95.45% of the total area and 


= !20< = accounts for 99.73% of the of the total area. This 


oO 
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implies that the filter can be truncated at 

w=t - to achieve shorter bandwidth with negligible error. 
€ 


For the region, ~2<w<®é the percentage error is 0.27% 
oO oO 


and the bandwidth is 24. 

Other properties of Gaussian filter in the frequency domain 
are shown in Figure 1. From the figure, it is observed that 
an increase in standard deviation leads to decrease in 
bandwidth of the filter. Figure 1 shows that standard 
deviation of 0.5 gives filter with the widest bandwidth 
while the standard deviation of 3 gives filter with the 
narrowest bandwidth. This implies that the bandwidth 
decreased with increase in standard deviation. On the other 
hand, smaller standard leads to higher amplitude. 
However, at the point where w=0, the amplitude of the filter 
is independent of the standard deviation. 

In the time domain, the Gaussian distribution is truncated 
at +30 since the amplitude of Gaussian kernel approaches 
zero when the execution time becomes more than 30. This 
can be shown in the illustrations in Table4. From Table 4, it 
is seen that at t = +40, the value of h(t) =0.0003 which is 
very close to zero. At any point when t = +5 0, the value of 
h(t) is zero. This implies that the filter can be truncated at 
t = +3 0 with very negligible error and shorter execution 
time. 

Although Gaussian filter can be truncated to enhance the 
speed of execution, the point of truncation must be 
carefully chosen to ensure minimum error. This error is 
calculated using the ratio of the normalized area excluded 
by —ko <t<ko region to the normalized area within 
-o<t<eo region. For a normalized 1-dimensional 
Gaussian filter the area within —oo < t < ~ region is unity 
(1). 

The percentage area occupied by a region determines the 
accuracy and execution time of the filter if it is truncated 
within the region. The areas within regions in Gaussian 
Table 1: Execution time of 3 x 1 Gaussian filter for 
different standard deviation computation expressions 


Standard deviation (o) | Execution time 
oa 2 
6 
gaa - | 080 3 
6 
_M+1) 067 4.04 
Sane | 


distribution are shown in Table 5. In Table 5, the area of a 
Gaussian distribution within the regions, -- < t <> and -50 
< t<5o are unity (1.0) each. 

The area occupied by the region - o < t < o accounts for 
68.27% of the total area of the distribution, the region -20 <t 
< 20 accounts for the 95.45% of the total area and -30 < t33. 0 
accounts for 99.73% of the of the total area. This implies that 
the filter can be truncated at -30 < ts 30 to achieve shorter 
execution time with a negligible error. For the region, -30 < 
ts 30 the percentage error is 0.27% and execution time is 
1.5. 

Other properties of the Gaussian filter in the time domain 
are shown in Figure 2. Figure 2 shows the plot of the 
amplitude of Gaussian filter in the time domain against 
time. From the figure, it is noticed that the filter with less 
standard deviation has the highest amplitude while a filter 
with higher standard deviation has the least amplitude at 
the point where t=0. Figure 2 also shows that filter with the 
highest amplitude take longer period to finish executing. 


1.5. Conclusion 


Based on the simulation carried out in this research, it is 
noticed that the standard deviation affects the performance 
and applications of Gaussian filter. The results show that a 
large standard deviation leads to long execution time and 
narrow frequency bandwidth of Gaussian filter. This 
implies that for application where speed is important, 
lower standard deviation is selected, but for application 
where a narrow frequency band is required for broad 
spectrum, the high standard deviation is required. It was 
observed that the error incurred in truncating Gaussian 
filter within a region in time and frequency domain is the 
same. However, the area enclosed by the same region is 


higher in the frequency domain. 


Table 2: Effects of frequency on the amplitude of 
Gaussian filter 


o =0.25 

H(0) =1 w= 0 
H(—4) = 0.6065 | H(4) = 0.6065 1 
w@=t- 
oO 
=p) = = 2 
H(-8) = 0.1353 | H(8) = 0.1353 Bee 
oO 
=12) = = 3 
H(-12) = 0.0111 | H(12) = 0.0111 —— 
oO 
=46) = = 4 
H(-16) = 0.0003 | H(16) = 0.0003 ease 
oO 
H(—20) = 0000 | H(20) = 0000 5 
wo=t- 
Je eG 
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Table 3: The error and bandwidth for regions of Gaussian 


Filter in frequency domain for o=0.25 
w Region | Region | Error (%) | Bandwidth Table 4: Effects of time on the amplitude of Gaussian 
Area Area (%) Frequency filter 
-o0<  < 00 10.0265 | 100 0.00 oo 
5 | 10.0265 | 100 0.00 40 il 
TZ SSG A(0) = 1.595 | t=0 
4 ey <4 | 10.0259 | 99.994 | 0.006 32 h(—1.25) = 0 h(1.25)=0 |t=+50 
(ome! h(—1) = 0.0003 h(1) = 0.0003 t=+40 
3. <3 | 99994 | 99.73 | 0.27 24 h(—0.75) = 0.0177 | h(0.75) = 0.0177 | t= +30 
nad, h(—0.5) = 0.2160 | h(0.5) = 0.2160 | t=+20 
ee gc Z | 95703 | 95.45 4.55 16 h(—0.25) = 0.9679 | h(0.25) = 0.9679 | t=+o 
er er 
te <i| 68450 | 68.27 31.73 8 
fer fer 
Table 5: The error and execution time for some regions in 
Gaussian filter 
t (o=0.25) Region | Region Error Execution 
Area Area (%) | (%) time 
-00< G < 00 100 0.00 °° 
5os<ts5o 100 0.00 2.5 
-4o <ts4o 0.99994 | 99.994 0.006 2.0 
30<t<30 | 0.9973 | 99.73 0.27 1.5 
2o<ts2o 0.9545 | 95.45 4.55 1.0 
-o<ts-o 0.6827 | 68.27 31.73 0.5 
——E—E——E——————————EE_ | 
~ Standard Devition=1 08 TX — 
» Standard Deviation=1.5 1? \ |—Standard Deviation=0.5 
08 + Standard Devation=2 0.7 7+ — Standard Deviation=1. 
a i creiveemlny al [|_| Standard Deviation=1.5 
£ eel $ = i | |-Standard Deviation=2 
206 05 | |__| — Standard Deviation=?.5 — 
3 3 f \--|-§tandard Deviation=3 
g £04 | 
20 Be 
; 4 0.3 
02 “y 
0.1 
a - 4 y 
Frequency (rad/sec) Time (sec) 
Figure 1: Plot of amplitude against frequency for -71< ws 7 Figure 2: Plot of amplitude against time for -7< ts 7 at 
at different standard deviation different standard deviation 
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